Effective gluon mass and 
infrared fixed point in QCD 



Arlene C. Aguilar and Joannis Papavassiliou 

Departamento de Fisica Teorica and IFIC, Centra Mixta, 
Universidad de Valencia - CSIC 
E-46100, Burjassot, Valencia, Spain 

Abstract. We report on a special type of solutions for the gluon propagator of pure QCD, 
obtained from the corresponding non-linear Schwinger-Dyson equation formulated in the 
Feynman gauge of the background field method. These solutions reach a finite value in the 
deep infrared and may be fitted using a massive propagator, with the crucial characteristic 
that the effective "mass" employed depends on the momentum transfer. Specifically, the 
gluon mass falls off as the inverse square of the momentum, as expected from the operator- 
product expansion. In addition, one may define a dimensionless quantity, which constitutes 
the generalization in a non-Abelian context of the universal QED effective charge. This 
strong effective charge displays asymptotic freedom in the ultraviolet whereas in the low- 
energy regime it freezes at a finite value, giving rise to an infrared fixed point for QCD. 
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A plethora of theoretical and phenomenological studies spanning more than two 
decades have corroborated the possibility of describing the infrared (IR) sector of 
QCD in terms of an effective gluon mass (for an extended list of references see 
According to this picture, even though the gluon is massless at the level of the 
fundamental Lagrangian, and remains massless to all order in perturbation theory, 
the non-perturbative QCD dynamics generate an effective, momentum-dependent 
mass, without affecting the local S'f/(3)c invariance, which remains intact j^]- 

The most standard way for studying such a non-perturbative effect in the 
continuum is the (appropriately truncated) Schwinger-Dyson equation (SDE) for 
the gluon propagator A^i,(g), defined (in the Feynman gauge) as 



P..(g)A(g^) + ^ 



9piu-'^- (1) 



Specifically, one looks for solutions having A(g^) reaching finite (no n- vanishing) 
values in the deep infrared, that may be fitted by "massive" propagators of the 
form A^^(g^) = + m?{q^). The crucial characteristic is that m?{q^) is not "hard", 
but depends non-trivially on the momentum transfer . When the renormalization- 
group logarithms are properly taken into account, one obtains in addition the non- 
perturbative generalization of g'^{q'^), the QCD running coupling (effective charge). 
The presence of m?{q^) in the argument of g'^{q^) tames the Landau singularity 
associated with the perturbative (3 function, and the resulting effective charge is 
asymptotically free in the ultraviolet (UV), "freezing" at a finite value in the IR. 



The running of m?'{q'^) is of central importance for the self-consistency of this 
approach, mainly because the value of A~^(0) is determined by integrals involving 
A(g^), m^(g^), and g'^{q^) over the entire range of (Euclidean) momenta. The UV 
convergence of these integrals depends crucially on how m^(g^) behaves as —>■ oo. 
If m?{q^) drops off asymptotically faster than a logarithm, then A~^(0) is finite. 
This, in turn, is crucial because the finiteness of A^^(O) guarantees essentially the 
renormalizability of QCD. 

In earlier studies of linear SDE [2I, [l| the m^i^q^) obtained drops in the deep 
UV as an inverse power of a logarithm. The main result reported in this talk is 
the existence of a new type of solutions for m?{q^) that drop asymptotically as an 
inverse power oi momentum (multiplied by logarithms) pi]. 



These solutions are found in the study of nonlinear SDE, in the framework 
defined from the combination of the Pinch Technique (PT) 0, 0, 0] and the 
Feynman gauge of the Background Field Method (BFM) [g, |7|, known as PT- 
BFM truncation scheme One of the most powerful features of the PT-BFM 
formalism is that, by virtue of the Abelian Ward identities satisfied by the various 
vertices, gluonic and ghost contributions are separately transverse, within each 
order in the "dressed- loop" expansion [ij]. This, in turn, allows one to truncate the 
series meaningfully, by considering only the diagrams (ai) and (a2) shown in FigH], 
(no ghosts included), without compromising the transversality of the answer. 

In order to reduce the algebraic complexity of the problem, we perform one 
additional approximation, dropping the longitudinal terms from the gluon propa- 
gators inside the integrals, i.e. we set A^p — *■ —igapA. Omitting these terms does 
not interfere with the transversality of the resulting propagator, provided that one 
drops, at the same time, the longitudinal pieces in the WI of Eq.dH) [l|,[3|. 

After these steps, the scalar function, A~^(g^) = g^ + in(g^), (where n(g^) is the 
gluon-self energy given by the diagrams (ai) and (a2) in Fig{T]) can be written as 



iP,,{q)A-\q^) = iV^,{q)q^ 



+ 




where the tree-level vertex F^q/j appearing in ([2]) is given by 



r^a/3(?,Pl,P2) = (pi-p2)M^a/3 + 2g/3^^„-2ga^M/3> 



(3) 



and ^va'p' represents the full three-gluon vertex. 
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FIGURE 1. The gluonic "one- loop dressed" contributions to the SDE. 



As a next step we will employ the "gauge technique" j^, expressing F as a 
functional of A, in such a way as to satisfy (by construction) the all-order Ward 
identity 

q>^f,^p{q,p,,P2) = z[A-l{p,)-A-'^{p2)], (4) 

characteristic of the PT-BFM. Specifically, we propose the following form for the 
vertex [3| 

r^"^ = Lf'^p + rr^ + , (5) 

with 



Tr^{q,PuP2) = -zc2{q^g^''-q"g^^) 



n(pi) , n(p: 



pI 



+ 



(6) 



Then, substituting Eqs.(l3]) and ([5]) into ([2]), introducing = x, = y, and 
defining the renormalization-group invariant quantity d{q'^) = g'^A{q'^) , we arrive 
at 



with 



i=l 



(7) 



Ai(x) 



As{x) 
Ai{x) 

Aq{x) 

As(x) 



l + ^]xl dyyC\y)d\y). 



Mix) = ^x / dyC{y)d{y), 



6c \ 

1 + -^] xC{x)d{x) J dyyC{y)d{y), 
+ ci^ C{x)d{x) dyy^C{y)d{y) , 

6c2 r 



5 Jo 



-C{x) 

5 X 



dyyC-{,y)d{y). 



dyy C{y)d{y). 



- I dyy'C\y)d\y). 
ox 



(8) 



The renormalization constant K' is fixed by the condition d ^(/i^) = /i^/fi'^, (with 
/i^ ^ A^), and /^(g^) = 61n(g^/A^), where A is QCD mass scale. Due to the poles 



contained in the Ansatz for F, d ^(0) does not vanish, and is given by the 
(divergent) expression 



57r2 



2(l + ci) d^kC{e)d{k')-{l + 2ci) / d^kk'C\k')d'ik') 



(9) 



which can be made finite using dimensional regularization, and assuming that 
m^^q^) drops sufficiently fast in the UV 

In order to determine the asymptotic behavior that Eq.Q predicts for m^{x) at 
large x, we perform the following replacements in the r.h.s. of ([S]) 

xC{x)d{x)^l, C{x)d{x)^l/x, C{y)d{y) = K{y), K{y) = \ (10) 

y + m^[y) 

Next, use the identity yA{y) = 1 —m'^{y)A{y) in all Ai{x), keeping only terms 
linear in . Then separate all contributions that go like x from those that go 
like on both sides, and match them up (oj. This gives rise to two independent 
equations, one for the "kinetic" term, which simply reproduces the asymptotic 
behavior xlnx on both sides, and an equation for the terms with m^(x), given by 



m^(x)lnx = C — ai dym^{y)K{y) -\ / dyyrn?{y)K{y) 

Jx ^ Jo 

+ ^ [ dyy^m\y)K{y) + a^x f dym\y)K\y), (11) 

^ Jo Jx 

with 

N 4 6ci 2 ^ 6c2 /^„s 

oi = -(I+C2-C1), 02 = - + — , 03 = --, 04 = 1 + — , (12) 

and 

_ /"OO 

C = b-^d-\0)+ai dym^{y)K{y). (13) 



Now, the important point to appreciate is that, in order for (ITTj) to have solutions 
vanishing in the UV, it is necessary to be sure that the constant term on the 
r.h.s. vanishes, i.e. C = 0. Since we know that (i~^(0) and the integral appearing 
in the r.h.s. of Eq.( |T3|) are manifestly positive quantities, it follows immediately 
that the C will be zero if and only if ai < 0. Notice that Eq. (fT3!) restricts the range 
of allowed values of the parameters ci and C2 through Eq. (|T2i) . In addition, and 
more importantly, it constrains the momentum dependence of m^(x) in the IR and 
intermediate regimes to be such that both terms on the r.h.s of Eq.f ll3l) cancel 
against each other. 

Assuming that Eq.f ll3l) is satisfied, it can be shown that Eq. fllll) admits the 
following asymptotic solutions for m {x) 0], 

ml{x) = \l{\nx)-^^+^'\ ml{x) = ^{\nxy^~\ (14) 



where Ai and A2 are two mass-scales, and 71 = — ai, 72 = 02- 

The first type of solutions, m?(a;), are familiar from studying linearized versions 
of Eq.([2]), see for example [2, llf- The second type of solutions, m^x), displaying 
power-law running, are particularly interesting because they are derived for the first 
time in the context of SDE. The possibility of an effective gluon mass dropping in 
the UV as an inverse power of the momentum was first conjectured in p|, and was 
subsequently obtained in the context of the operator-product expansion [lO| ; there 
the resulting gluon self-energy was identified as the effective gluon mass, leading 
to the relation m^(a;) ~ {G^)/x, where (G^) is the dimension four gauge-invariant 
gluon condensate. 

Which of the two types of solution will be actually realized depends on the details 
of the three-gluon vertex, F, and more specifically on the values of the parameters 
Ci and C2. Our numerical analysis reveals that the sets of values for ci and C2 giving 
rise to logarithmic running belong to an interval that is disjoint and well-separated 
from that producing power-law running. In what follows we will focus our attention 
on the latter type of solutions. In Figj2]we present typical solutions for the (i(g^), 
m^(g^) and the effective charge a(g^) = g^{q'^)/AiT. 
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FIGURE 2. Upper panel: the numerical solution for d{(f'). Lower panels: Left: dynam- 
ical mass with power-law running, for = 0.5 GeV^ and p = 1.046 in Eq. (jl7p . Right: 
the running charge, a{q^) = g^{q^)/A'K, fitted by Eqs. (fT5|l and ([TBI) . 



The way to extract from d{q^) the corresponding m?{q^) and g'^{q^) is by casting 
the numerical solutions shown in Fig|2] into the form 



q2 _|_ jy^2(^q2\ 



bin 



2 ^2^^2\\ \ 1 -1 



q^ + f{q^,m ^{q^)) 
A2 



with 



' q^ + m?[q'^) [q^ + m? [q^)]"^ 

where pi, p2, and ps are fitting parameters. 

The functional form used for the running mass is 



(15) 



(16) 



2/ 2\ 

m [q 



■nin 



q^ + 



A2 



A2 



72-1 



(17) 



In the deep UV Eq. ffTTl) goes over to mKg^), whereas at = it reaches the finite 



value m (0) 



Note that f {(f ,rr? {(f)) is such that /(0,m^(0)) > 0; as a result 



the perturbative Landau pole in the running coupling is tamed, and q^fq ^) reaches 
a finite positive value at = 0, leading to an infrared fixed point [2I.T11I. 1^ . 

To summarize our results, from a gauge-invariant SDE for the gluon propagator 
we have derived an integral equation that describes the running of the effective 
gluon mass in the UV, and have demonstrated that, depending on the values of 
two basic parameters appearing in the three-gluon vertex, one finds solutions that 
drop as inverse powers of a logarithm of g^, or much faster, as an inverse power of 
g2. Moreover, we have extracted an asymptotically free effective (running) charge, 
that freezes in the low-momentum region, implying the existence of a IR fixed point 
for QCD. 
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